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The ground-state properties of a spin 5 = 1/2 tetrameric Heisenberg antiferromagnetic chain 
with alternating couplings AF1-AF2-AF1-F (AF and F denote antiferromagnetic and ferromagnetic 
couplings, respectively) are studied by means of the density matrix renormalization group method. 
Two plateaux of magnetization m are found at m = and 1/4. It is shown that in such a spin-1/2 
AF system, there is a gap from the singlet ground state to the triplet excited states in the absence 
of a magnetic field. The spin-spin correlation function decays exponentially, and the gapped states 
have a nonvanishing string order, which measures a hidden symmetry in the system. By a dual 
transformation, the string order is transformed into a ferromagnetic order and the hidden symmetry 
is unveiled to be a Z2 x Z2 discrete symmetry, which is fully broken in the gapped states. This 
half-integer spin Heisenberg AF chain is in a Haldane-like phase, suggesting that the present findings 
extend the substance of Haldane's scenario. A valence-bond-solid state picture is also proposed for 
the gapped states. 
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I. INTRODUCTION 

Even though many models have been extensively stud- 
ied experimentally and theoretically, one-dimensional 
quantum spin systems are still an attractive field in low- 
dimensional quantum magnetism. One of motivations is 
from Haldane's conjecture (l|. Haldane proposed that 
an isotropic Heisenberg antiferromagnetic chain (HAFC) 
with integer spin has a finite gap from the singlet ground 
state to the triplet excited states, and the spin-spin cor- 
relation function decays exponentially; while the HAFC 
with half-integer spin has a gapless spectrum and a cor- 
relation function with power-law decay. Although Hal- 
dane's conjecture has been confirmed experimentally and 
numerically in many systems, there is not a rigorous 
proof till now. By means of the valence bond states, 
a rigorous disordered ground state for the biquadratic 
Heisenberg Hamiltonian with spin 5 = 1 was proposed 
by Affleck, Kennedy, Lieb and Tasaki (AKLT) ||. This 
AKLT model was shown to have a spin gap from the 
singlet ground state to the triplet excited states, an ex- 
ponentially decaying correlation function, and a nonlocal 
topological order which is different from the long-range 
order of spin-spin correlation, thus confirming the Hal- 
dane's conjecture based on a specific model. This non- 
local string order, which is regarded as a hidden anti- 
ferromagnetic Neel order, was found by Den Nijs and 
Rommelse Q in the spin-1 HAFC. Kennedy and Tasaki 
[3| introduced a nonlocal unitary transformation to re- 
veal the hidden symmetry for this string order, which 
was found to be a discrete Z2 x Z2 symmetry. The Hal- 
dane phase of integer-spin HAFC is thus characterized 
by the complete breaking of the Z2 x Zi symmetry Q ■ 

Besides the isotropic HAFC with integer spin, Haldane 
phase has also been found in other spin systems, like some 
5 = 1/2 spin ladders H,H,0| and spin-1/2 ferromagnetic- 
antiferromagnetic (F-AF) alternating Heisenberg chain 



Q. In these systems, the gap and the string order vary 
monotonically with the ratio between F and AF interac- 
tions {Jf/Jaf) and when Jp — > —00, recover the val- 
ues obtained for the spin-1 HAFC. The weak F coupling 
phase belongs to the same phase as the strong F coupling 
case. When the F couplings dominate the AF ones, the 
dimers of two spins coupled by F interactions behave as 
spin-1, and the systems reduce to the spin-1 HAFC. The 
Z2 x Z2 symmetry is also fully broken in these systems, 
which indicates that the systems belong to the Haldane 
phase. 

In this paper, we are concerned about the existence 
of a Haldane-like phase in an AF chain with half-integer 
spin, which, unlike the spin-1/2 F-AF alternating chain, 
could not be reduced to a HAFC with integer spin in any 
circumstance. The spin-1/2 trimerized F-F-AF Heisen- 
berg chain might be a choice, but its spectrum is gapless 
Q. Here we consider a spin 5=1/2 tetrameric HAFC 
with alternating couplings AF1-AF2-AF1-F. As far as we 
know, there has been no report on any tetrameric an- 
tiferromagnet yet, and only has a tetrameric ferrimag- 
netic Heisenberg chain compound Cv^-ClpyWNj^ (F- 
F-AF-AF) been widely studied ES M, Gj, El El El • 
By means of the density matrix renormalization group 
(DMRG) method and dual transformation, we have 
found that this spin-1/2 tetrameric system is in a gapped 
phase, most properties of which are compared with the 
features of the Haldane phase. Most importantly, the 
phase we found has a string order and a hidden Z2 x Z2 
symmetry, which is also fully broken in this system. 
This alternating tetrameric HAFC with spin-1/2 is in a 
Haldane-like gapped phase. In this sense, our findings ex- 
tend the substance of Haldane's scenario, which implies 
that the spin gap can also appear in certain half-integer 
spin HAFCs. 

The rest of this paper is organized as follows. In Sec. 
II, we shall introduce the model Hamiltonian, and give 
a brief discussion on the isolated tetramer systems. In 
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FIG. 1: (Color online) Schematic representation of the 
tetrameric Heisenberg spin chain with alternating couplings 
described by Hamiltonian Jl]). 

Sec. Ill, we shall present our DMRG results on the mag- 
netic properties of the system in a longitudinal magnetic 
field, as well as the gap behaviors, spin-spin correlation 
function and string order in the gapped ground states in 
absence of the external field. A dual transformation is 
introduced to the model in Sec. IV to unveil the hidden 
symmetry. In Sec. V, we shall propose a valence-bond- 
solid state picture and a trial wavefunction for the gapped 
states. Finally, a summary and discussion will be given 
in Sec. VI. 



II. SPIN-1/2 TETRAMERIC HEISENBERG 
ANTIFERROMAGNETIC CHAIN 

The Hamiltonian of the isotropic tetrameric HAFC 
with alternating couplings AF1-AF2-AF1-F in a longi- 
tudinal magnetic field is given by 

H = y^(J">tjr 1 S4j_3 • S4./-2 + ■ S4J-1 

3 

+ JaF^^-I ■ S4,' — JF$4j ■ S4J+1) — h/lSj,(l) 

3 

where Jaf x 2 (> 0) denote the different AF couplings, 
Jf (> 0) denotes the F coupling, and h is the external 
magnetic field. We take gfiB — 1 for convenience. The 
schematic spin configuration for the tetrameric HAFC is 
depicted in Fig. [1] Here we would like to remark that 
a spin-5/2 tetrameric HAFC compound, C44H 36 Ni 6 Mn 2 , 
with alternating couplings AF1-AF2-AF1-F was synthe- 
sized experimentally [l6j]. Since the DMRG calculations 
on most physical properties of such a spin-5 /2 tetrameric 
HAFC are impermissible for our present computing ca- 
pacity, and considering that the spin-1/2 tetrameric 
HAFC is more fundamental and interesting, and can be 
readily made comparisons to other Haldane spin systems, 
we opt to pay attention to the tetrameric system with 
spin-1/2 in the following. Meanwhile, for the reasons 
of comparison, some results on the spin-5/2 tetrameric 
HAFC will also be calculated. 

Let us first discuss some cases of the present model. 
When JaFx = 0, the system is decoupled into isolated 
dimers with Jaf 2 or Jf couplings. The dimers coupled 
by J af 2 form local singlets while those coupled by Jf 



form local triplets. Thus, the excitation of the system 
from S* ot = to S* ot — ±1 is gapless due to the triplets. 

When Jaf 2 — 0; the system is decoupled into isolated 
tetramers, which becomes exactly soluble. In the absence 
of an external field, the ground state of the tetramer is 
antiferromagnetic with Stot = 0, and the first excited 
state is threefold degenerate with Stot — 1 • The gap A is 

A = \{^J\ Fl + 2J AFi Jf + J 2 f~ Jf), (2) 

which increases monotonically with J af^ but decreases 
with Jp ■ When Jp 3> Jaf 1 , the gap has an asymptotic 
form A ~ JaFxQ + 3Ja_Fi/2J^ )/2 and equals J^Fi/2 
when J f — * 00 . 

In the limit Jf = 0, the ground state of the tetramer is 
also antiferromagnetic with Stot — 0, and the first excited 
state is threefold degenerate with Stot = L The gap A 
from the ground state to the first excited state is 



A = 2^ JaFi + V 4J ^ Fl ~ 2J afiJaf 2 + J\f 2 

~ \J J\f x + J\f 2 )i ( 3 ) 

which enhances monotonically with J af\ but decreases 
with Jaf 2 ■ When Jaf 2 Jaf x , the gap behaves as 
A ~ J\ F JU AF-i and vanishes when Jaf 2 — * 00. 

In the large Jp limit, the two spins coupled by Jp 
form a spin-1, and the model reduces to a Heisenberg 
chain with alternating spin-(l,l/2,l/2): 

Hj F ^oo — y^(^JAFi^3j-2 ' S3j_i + JAF 2 ^3j-l ' &3j 

3 

+ lj AFi S 3r S 3]+1 ), (4) 

where Sj is the spin operator with S = 1. This model 
appears to be rarely studied. 

For arbitrary couplings, we will apply the DMRG 
[l7l [lli ] method to explore the ground-state properties 
of the tetrameric chain defined in Eq. ((T|). In the fol- 
lowing DMRG calculations, the chain length is taken as 
N = 160, and the Hilbert space is truncated to 100 most 
relevant states. In the calculations for correlation func- 
tion and string order, 160 optimal states are kept for 
accuracy. Open boundary conditions are adopted. The 
truncation error is less than 10 -8 in all calculations. 



III. MAGNETIZATION, HALDANE-LIKE GAP, 
CORRELATION FUNCTION AND STRING 
ORDER 

A. Magnetization 

Due to the competition between interactions and the 
external magnetic field, a quantum magnet often shows 
exotic properties under magnetic fields. A field can 
close the zero-field gap, and may induce magnetization 
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plateaux under some conditions [ljj [2(| [2l|, ■ There- 
fore, we shall investigate the magnetic properties of the 
present alternating tetrameric system in a longitudinal 
field. Many cases with different couplings have been cal- 
culated. Here the results of the case with Jaf x '■ Jaf 2 '■ 
Jf = 1 : 1 : 1 are presented for example. Fig. [2] (a) 
shows the magnetic field dependence of the magnetiza- 
tion per site m(h). The system exhibits two magnetiza- 
tion plateaux, a zero-field gap and a m = 1/4 plateau, 
both of which satisfy the Oshikawa-Yamanaka- Affleck 
(OYA) [2lJ condition n(S — to) = integer, with n the 
period of the ground state, S the magnitude of spin and 
to the magnetization per site. A similar magnetization 
curve has been found in the p-merized antiferromagnetic 
spin chain with p = 4 [23[ . Both systems have shown all 
plateaux that the OYA condition permits. The to = 
plateau in the present system exhibits a spin gap in the 
spectrum, which is from the singlet ground state to the 
triplet excited states. Such a gap has also been found in 
the spin- 1/2 F-AF alternating Heisenberg chain Q. In 
such spin-1/2 chains with even lattice translation sym- 
metry, the spectrum could have a gap even if the ground 
state has no spontaneous translation symmetry breaking 
[2ll |24| . In contrast, spin-1/2 chains with odd lattice 
sites per unit cell are gapless if the translation symmetry 
does not break spontaneously, like the spin-1/2 Heisen- 
berg chain [25| and the spin-1/2 F-F-AF trimerized chain 
. The further discussions on this gapped ground states 
will be presented below. 

For the to = 1/4 plateau, the local magnetic moment 
(Sj) and the spin-spin correlation function (SqS*) are 
calculated. These two quantities show a perfect sequence 
with a period of 4 and (S^S*) has a long-range order. It 
is shown that (Sj) behaves as {• • • , (Si, S%, S 2 , Si), ■ ■ ■} 
with S x = 0.4468 and S 2 = 0.0532, giving rise to the 
magnetization per site m = 1/4, as shown in Fig. [2] (b). 
The plateau state can be described by an approximate 
wave function, 

Ipi = a \ T4i-3t4i-2t4j-ll4i) + b\ t-lj-3 ?4i — 2j-4i— 1 ?4i) 
+ c\ T4i-3l4i-2t4j-lt4i) + d\ J.4i_3 T4i-2?4i- 1 ?4i) 

(i = l,---,N/A), (5) 

where fj (lj) denotes spin up (down) on site j and N is 
the total site number. If the coefficients a 2 = d 2 = S2 and 
b 2 = c 2 = S\, the local magnetic moment and the spin- 
spin correlation function deduced from this wave function 
perfectly fit into the numerical results. The formation 
of the to = 1/4 plateau in the limit Jafi = is quite 
simple, which is helpful for understanding the general 
case. When JyiFi = 0, the system is decoupled into the 
isolated dimers coupled by Jaf 2 or Jf- After applying 
a magnetic field, the dimers coupled by Jf, which are 
in the triplet states, become polarized and, the system 
exhibits a magnetization plateau at m = 1/4. The width 
of the plateau is determined by the upper critical field 
J AFi leading to full polarization. 

The magnetization process at h Cl < h < h c . 2 and h C3 < 
h < h s are shown in Figs. [U(c) and (d), respectively. It is 
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FIG. 2: (Color online) (a) The external field dependence of 
the magnetization per site m(h) for the spin-1/2 tetrameric 
HAFC with couplings Jaf 1 :Jaf 2 'Jf=^'^ (b) The local 
magnetic moment (5*1) in the magnetization plateau states. 
The DMRG results of m(h) in the gapless states can be fairly 
fitted by Eq. (|6]) for (c) h ci < h < h C2 , and (d) h C3 < h < h s . 

found that the behaviors of the curves could be described 
as 

m(h) = mi + (h - hi)(m 2 - mi)/(h 2 - hi) 

+ a(h 2 -h)y/h-hi- (3(h - h 1 )yjh 2 - h, (6) 

where hi (h 2 ) and mi (m 2 ) are the lower (upper) critical 
field and the corresponding magnetization per site, a and 
(3 are two parameters. For h Cl < h < h C2 , a = 0.95 and 
(3 = 1.2, while for h C3 < h < h s , a = 3.2 and (3 = 4.0. 
The curves obtained from Eq. ([6|) are also shown in the 
figures, which fit into the DMRG results quite well. 

B. Haldane-like gap 

To understand the properties of the ground state of 
the spin-1/2 tetrameric alternating HAFC in the absence 
of a magnetic field, the spin gap, spin-spin correlation 
function, and string order are studied by means of the 
DMRG method in the whole parameter region, which 
will be compared with the features of Haldane phase. 

We consider the spin gap A from the ground state to 
the triplet excited states, namely, 

A = Ei - E , (7) 

where Eq is the ground-state energy and Ei is the low- 
est energy in the subspace with Sf ot = 1, and Jaf x is 
chosen as the energy scale. Jf and Jaf 2 dependences 
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of the gap are presented in Figs. [3J (a) and (b), respec- 
tively. Fig. [3] (a) shows the gap as a function of Jp for 
Jaf 2 = 0.0, 0.1, 0.5, 1.0 and 2.5. When J A f 2 = 0, the gap 
is determined by Eq. (J2J). It is shown that with the in- 
crease of Jf, the gap A smoothly decreases and converges 
to a certain value A c when Jf — > oo. A c is determined 
by Jaf 2 - I n the limit of Jaf 2 — 0, A c = J^p 1 /2. In 
the large Jf limit, the system is equivalent to the spin- 
(1,1/2,1/2) model, which will be proved to have a gap 
by the non-linear a model below, where A c is related 
to the gap of this model. With the increase of J F , the 
tetrameric system changes from a tetramer model to the 
spin-(l, 1/2,1/2) model, and the gap decreases continu- 
ously. In Fig. [3] (b) , the gap is plotted as a function of 
Jaf 2 for J F = 0.0,0.1,0.5,1.0 and 2.5. When J F = 0, 
the gap is evaluated by Eq. ([3]). It is observed that the 
gap decreases rapidly with increasing J af 2 and tends to 
zero when Jaf 2 — ► oo. With the increase of Jaf 2 , the 
tetrameric chain changes from a tetramer model to a sys- 
tem with relative small Jaf x interactions, which is gap- 
less when Jaf ± = 0. From the above discussions, it can 
be seen that both Jp and J af 2 give rise to the decrease 
of the gap, but with different behaviors. The influence 
of JaFx on the gap is found to be distinct from those of 
J af 2 and Jp ■ The gap would increase with Jaf x ■ Be- 
cause the gap only vanishes when Jaf 2 / 1 JaFx — > oo, it is 
reasonable to conclude that a gap would be generated by 
an arbitrary small Jaf x ■ 

Next, let us consider the asymptotic behaviors of the 
gap. In Fig. [3J (c), the gap as a function of JaFx/Jf 
is plotted in large Jf limit for J af 2 =0.5 and 1.0. It 
is shown that the gaps linearly converge to A c . J af 2 
alters the value of A c but does not change this linear 
asymptotic behavior. Thus, the variations of the gap for 
large Jp could be evaluated by 



A = J AFl ( 



JAFx 



JAF! x 



(8) 



where A c and the parameter fi are determined by Jaf 2 - 
In the large J af 2 limit, the isolated tetramer system with 
Jp = gives A ~ J\ f JJaf 2 - The inset of Fig. [3J (c) 
shows the gap behaviors for Jp — 1.0 and 2.5 for large 
J af 2 i which exhibits an asymptotic behavior 
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A r~~l AFl C ~* J AF, /JaFo 

Jaf 2 



(9) 



where 7 is a parameter, and 7 = 3.2 for Jp = 1.0 and 
1.8 for Jp = 2.5. The extrapolations show the linear 
behaviors near Jaf 1 — 0. When J A f 1 dominates, the 
gap behaviors are shown in Fig. [3J (d). It can be seen 
that the gap decreases linearly with increasing both J af 2 
and Jf for large Jaf ± ■ 

The DMRG results suggest that the system is in disor- 
dered spin liquid states with a spin gap. The gap varies 
monotonically with the couplings, and does not show any 
singularity. We have calculated the ground state energy 
against the couplings and found that the ground state en- 
ergy varies also monotonically with the interactions, and 
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FIG. 3: (Color online) (a) The Jp dependence of the gap A 
for Jaf 2 =0.0, 0.1, 0.5, 1.0 and 2.5. (b) The Jaf 2 dependence 
of the gap A for ^=0.0,0.1,0.5,1.0 and 2.5. (c) The gap 
behaviors with Jaf 1 /Jf for large Jf- The inset shows the 
gap behaviors with Jaf 1 /Jaf 2 for large Jaf 2 - (d) The gap 
behaviors with Jf for large Jaf x ■ The inset shows the gap 
behaviors with Jaf 2 for large Jaf ± ■ 



there is no any nonanalyticity. This nonanalyticity of the 
ground state energy suggests an absence of the quantum 
phase transition in this system [26| . 

It has been shown that the gap always exists if JaFx 
is finite. However, it is hard to identify the valid- 
ity of this statement from numerical calculations when 
Jaf 2 3> J AFi- Thus, we consider the nonlinear a model 
(NLSM) of this tetrameric system. Because the gap de- 
creases with Jp, we consider the case in the limit of 
Jf — > 00, when this tetrameric system is equivalent to 
the spin-(l,l/2,l/2) model which is described by Eq. 
The NLSM for such spin chains has been develope d by 
Affleck [13, Fukui and Kawakami [H, and Takano 
The magnitude of Sj is denoted as Sj, which could be 1/2 
or 1 here. The period of the system 2b is regarded as 6 
for convenience. After a standard procedure [29(, we have 
the effective action for the spin-(l,l/2,l/2) Hamiltonian: 
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eff 
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pL 

dT J dx{ 2ajW jW J(D 

.J[ 0] ,„ „ s ajW 
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— -m • (<9 T m x d x m) 



(^m)^},(10) 



where 1//3 denotes the temperature, a is the lattice spac- 
ing, L is the length of the chain, m presents the spin 
variables, and {J™} are defined by 
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FIG. 4: (Color online) The gap A of the spin-5/2 tetrameric 
HAFC as a function of the inverse of chain length (1/L) for 
(a) J F = 1.0, J A f 2 =1.0,2.0 and 5.0; and (b) J AF2 = 1.0, 
J F = 1.0, 2.0 and 5.0. 



with J a 



Jo 



and 



= ELi(-lf +1 ^ The 
topological angle 8 is 47T J^ / and the NLSM provides 
a gapless condition when Q/2ir is a half-odd integer. After 
applying these equations to the spin-(l,l/2,l/2) model, 
we have the gapless condition 



J AFl 61 - 11 



Jai 



6 — 4Z 



(I = positive integer). 



(12) 



As Jafx and Ja_f 2 are both positive, Eq. (fT2")l could not 
be satisfied in any case, which indicates that even in the 
large Jp limit, the present tetrameric chain is gapped 
for any Jaf x > and Jaf 2 > 0. Considering the ef- 
fect of Jf on the gap, we could conclude that this spin- 
1/2 tetrameric Heisenberg AF chain is always gapped if 
J AFl £ 0. 

For a comparison, the gap of the spin-5/2 tetrameric 
HAFC is computed by utilizing the DMRG method. The 
chain length is taken as 140 and the Hilbert space is trun- 
cated to 400 optimal states. The truncation error is less 
than 10 -13 in energy calculations. In Fig. 0] the DMRG 
results show that the extrapolations of the gap converge 
to nonzero values. Similar to the case with spin- 1/2, the 
gap of this spin-5/2 tetrameric HAFC also decreases with 
both increasing Jaf 2 an d Jf- The gap diminishes more 
rapidly with increasing Jaf 2 > which is also a feature of 
the spin-1/2 case. Although the gap of the spin-5/2 sys- 
tem is smaller than that in the spin-1/2 case, the qual- 
itative behaviors of the gap for both systems appear to 
be similar. 



C. Correlation function and static structure factor 

To characterize the gapped ground states of the 
tetrameric spin chain, let us investigate the spin-spin cor- 
relation function (SfS?) and the static structure factor 
S(q) that is defined as the Fourier transform of the cor- 
relation function: 
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FIG. 5: (Color online) The DMRG results of the static struc- 
ture factor S(q) in the gapped ground states of the spin-1/2 
tetrameric HAFC for (a) J F =0.1, J AF2 =0.1, 1.0 and 5.0; (b) 
Jf=5.0, Jaf 2 =0.1, 1.0 and 5.0; (c) J F =0.0, Jaf 2 =0.1, 1.0 and 
5.0; (d) Jaf 2 =0.0, J f =0.1, 1.0 and 5.0. 



Our DMRG results show that the correlation function 
decays exponentially in the gapped states. For different 
couplings, the correlation length and the behavior of the 
correlation function have dramatic changes, which could 
be displayed more clearly in the static structure factor 
S(q). 

We have performed a large amount of calculations on 
the correlation function and static structure factor in a 
wide range of the parameter region. The changes of S(q) 
are shown in Figs. [5] (a) and (b) for some parameters as 
examples. In Fig. [5] (a), for Jp =0.1, Jaf 2 = 0.1, 1.0 and 
5.0, S(q) shows an obvious maximum at q — it, and with 
the enhancement of Jaf 2 , there are two small maxima 
appearing at q = ir/3 and 5ir/3. When Jaf 2 — 0.1 and 
Jf = 5.0, the maximum at q — tt becomes a valley and 
two small maxima appear near q = 37r/5 and 7ir/5. With 
further increasing Jaf 2 , S(q) exhibits a totally different 
behavior, showing four peaks at q = tt/4, 37r/4, 57r/4 and 
77r/4, as shown in Fig. [5](b). In these numerical results, 
S(q) shows three different behaviors due to the competi- 
tion of the couplings. It will be shown that some features 
of S(q) mentioned above are determined by the short- 
range correlations, and some others are due to the in- 
crease of the correlation length and the translation sym- 
metry of the ground states. 

The correlation function of this model has no long- 
range order, but exhibits a short-range order charac- 
terized by some q m ax of S(q). In the limit Jaf 2 = 
and Jf = 0, the tetrameric chain is decoupled into sin- 
glet dimers, and the static structure factor has a sim- 
ple form S(q) = (1 — cosq)/4. The short-range correla- 
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tions (S'fj-i'Sfi) are reflected by the maximum of S(q) 
at q — ir. When Jaf 2 is switched on, the tetrameric 
chain becomes a tetramer system. The short range cor- 
relations of the tetramer produce another two maxima of 
S(q) at q = ir/3 and 5ir/3, as shown in Fig. \5\ (c). On 
the other hand, when Jf is turned on, the short-range 
correlations drive the maximum of S(q) at q = n to split 
into two small maxima, which move towards g=37r/5 and 
77r/5 with increasing Jp, as presented in Fig. [5] (d). The 
distinct behaviors of S(q) in the two situations are due 
to the different couplings of the tetramers. These results 
found in the tetramer systems also explain the behaviors 
of S(q) for small Jaf 2 ° r Jf shown in Figs. [5] (a) and 
(b). When Jaf 2 or Jf is smaller than Jaf ± , the correla- 
tion length of the system is quite short, reflecting a highly 
disordered ground state, and the correlation function and 
static structure factor could be well characterized by the 
corresponding decoupled tetramer system. 

When both Jaf 2 and Jf are large enough, S(q) ex- 
hibits a totally different behavior. There are four peaks 
at q = 7r/4, 37r/4, 5tt/4 and 77r/4, as shown in Fig. [5] 
(b) when Jaf 2 = Jf = 5.0. In this case, the correlation 
function keeps exponentially decaying but shows a struc- 
ture with a period of 4, and its correlation length becomes 
large. It is also observed that the correlation function has 
a characteristic such as {SfSj) — {S*S* +1 } if the spins Sj 
and Sj+i are coupled by J F , and (SfS z ) ~ — (5f 5J +1 ) 
if the spins Sj and Sj+i are coupled by Jaf 2 ■ This fea- 
ture reflects that the spins coupled by Jf behave like 
triplets, and the spins coupled by J af 2 behave like sin- 
glets. As the correlation function that decays exponen- 
tially is modulated with a period of 4, the static struc- 
ture factor S(q) shows the peaks at q = ir/4, 37r/4, 57r/4 
and 77r/4. Such a modulated structure of S(q) has 
also been observed in the spin-1/2 alternating F-AF 1 30] 
and spin-1/2 trimerized F-F-AF Heisenberg chains [3l| . 
which show the peaks of S(q) at q = ir/2,3ir/2 and 
q = 7r/3, 7r, 57t/3, respectively. Generally speaking, for an 
exponentially decaying correlation function that is mod- 
ulated with a period of n, the static structure factor S(q) 
could show antiferromagnetic peaks at q = tt(1 + 21) /n 
with I = 0, 1,2, • • • ,n - 1. 

Although the ground state energy has no nonanalyt- 
icity in the parameter space, namely, no quantum phase 
transition happens, the static structure factor S(q) shows 
rich features of the ground states for different couplings 
in this tetrameric system. 



in position, their sequence has a hidden Neel order, that 
is, if all sites with Sf — are removed, the left sites 
with Sf — 1,-1 have a Neel order. The string order 
was found to be a common feature of Haldane phase. 
Hida suggested a string order for spin-1/2 F-AF chain 
[1] , which was used to characterize the phase diagram of 
this spin chain [32| . The string order for spin ladders has 
also been considered in distinct gapped phases [33j. 

Here we suggest a string order for this tetrameric sys- 
tem, and use it to characterize the disordered spin liquid 
phase. Considering the basic feature of string order and 
the translation symmetry of the tetrameric Hamiltonian, 
we define four string operators for this tetrameric chain, 
which are denoted by 6£_ 3i4i+2 , 0£_ 2i4i+3 , 
and ®4i-3.4j+2 i s defined as 



^41-3,4^+2 



-SS_ 3 exp(*7r V .S'l.S' 



fc=4i-2 



a \ net 
k )°4j+2i 



(14) 



and the corresponding order parameter is 



^4i-3,4i+2 — . . U P (©4i-3.4i+2) i (15) 
|j-t|-»oo 

where i and j denote the unit cell and a = x,y ot z. 
The other three ones are defined in the similar way and 



they could be obtained by a translation of 0?, 



4i-3,4j+2- 



As 



the present model in the absence of a magnetic field has 
a SU(2) symmetry, we need only to consider Q z . The 
spin configuration shows that Ofj_ 34j+2 and 0|^_i4j 



have the same topological structure, while 0|j_ 



4i-2,4j+3 



and 



4i 4 - +1 possess another same structure. The DMRG cal- 
culations have been performed on the four string orders. 
It was found that the string orders of the first structure 
have the same finite value, while the string orders of the 
second structure are both zero in the whole parameter 
space. For the spin-1/2 F-AF alternating chain, the two 
different definitions of string order also have two different 
values, one finite and one zero (34[. This difference would 
be explained in the next section. Here we only consider 
the nonvanishing string order, which would be denoted 
as O z strmg below. 

When both J af 2 and Jf are equal to zero, the string 
order has its maximum 0.25, as shown in Fig. [BJ Af- 
ter tuning up Jf or Jaf 2 , the system becomes tetramer 
assemblies. In the limit J af 2 = 0, the string order is 
evaluated from the wavefunction of the tetramer ground 
state as a function of Jp / Jaf^ '■ 



D. Topological string order Q 



string 



In order to investigate natures of the Haldane-like 
phase in the disordered phase of the tetrameric chain, 
we have calculated the string order to detect the hid- 
den symmetry in the gapped phase. The nonlocal string 
order was found in spin-1 antiferromagnetic Heisenberg 
chain by Den Nijs and Rommelse They found that O z string 
although the spins with Sf = 1,0,-1 are not ordered 



4+ 7^ + ^4 + 2^- + (^ 
12 J 4 + 2-^- + (^ 



(16) 



while for Jp = 0, the string order as a function of 
J af 2 I Jaf! is expressed as 



^ _ Jaf 2 
Jaf± 



4-2 



Jaf 



F 2 I I J AF 2 \ 2 
AFx ^ JAF 1 ^ 



^/->-2^ + (^) 2 



(17) 
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as shown with solid lines in Figs. [6] (a) and (b), respec- 
tively. From these equations and numerical simulations, 
one may find that the string order in the case of Jaf 2 = 
has an asymptotic behavior: 



O z J-fJ^l 2 

^string a a a V j ) 



64 "J AFl 

when Jf -C Jaf x ■ In the case of Jf — 0, it behaves as 



(18) 



n z 1 1 , JaF 2 

string ~ 4 64 



(19) 



when J^4f 2 ^ >7a Fi ■ On the other hand, in the limit of 
Jaf 2 — 0, the string order has an asymptotic form when 
Jf ^> Jaf x , 



O 



string 



while it behaves as 



1 J AFx 

4 J F 



"2 J r AF 1 i Jf 



O 



string 



1 Jaf x 
4 Jaf 2 ' 



(20) 



(21) 



when Jaf 2 3> JaFi for the case Jp = 0. 

Besides the limiting cases, the string order is evaluated 
as a function of the couplings by means of the DMRG 
method to detect the hidden symmetry in the gapped 
phase. Fig. [6] presents the DMRG results of the string 
order as functions of Jf/Jaf 1 and Jaf 2 / Jaf x - In Fig. [6] 
(a), the string order is shown as a function of Jf / ' Jaf 1 for 
Jaf 2 = 0.1,0.5,1.0 and 2.5. The string order decreases 
monotonically with Jp/ Jaf ± and converges to a nonzero 
value when Jf/Jafx — > oo. In Fig. [5] (b), the string 
order as a function of Jaf 2 / 'Jafi is evaluated for Jp = 
0.1,0.5, 1.0 and 2.5. It is shown that it decreases rapidly 
with J af 2 I JaFi and vanishes when Jaf 2 /JaFi —* oo. 
As is seen, the string order is always finite in the gapped 
phase, indicating the existence of a hidden symmetry in 
this phase. The variations of the string order show the 
same tendencies as the gap behaviors. The string order 
has its maximum when Jp — and J af 2 = 0, where the 
system also shows the maximum of the gap. When the 
gap vanishes in the case of Jaf x = 0, the string order is 
also zero. This indicates that the opening of the gap and 
its behaviors could be characterized by this string order. 

The asymptotic behaviors of the string order are shown 
in Figs. [6] (c) and (d). Fig. [6] (c) shows the string orders 
when Jf^$>Jaf! for J af 2 =0.5 and 1.0. The inset shows 
the string orders when Jaf 2 ^>JaFi for Jf — 0.5 and 1.0. 
Fig. [5] (d) and the inset show the string orders when 
Jaf^Jf and Jaf 1 ^ > Jaf 2 , respectively. It is found that 
the asymptotic behaviors of the string orders preserve the 
features as those of the corresponding tetramer models, 
which are described by Eqs. (fI5 ]) -(|2"T| ) . 

It should be noticed that although the string order is 
nonlocal, it could measure some localized singlet corre- 
lation that depicts the singlet state. The string order 
has been found to characterize the behaviors of the gap 
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FIG. 6: (Color online) The string order as a function of 
Jf/Jafx (a) and Jaf 2 / Jaf x (b). The asymptotic behavior 
of the string order for large Jf or Jaf 2 (c) and for large Jaf x 
(d). 



perfectly, as shown in Figs. El (a) and (b). The gap mono- 
tonically increases with Jafi > and is related to the singlet 
states of the spins coupled by Jaf x , which could also be 
measured by the string order. When the system is the 
assembly of localized singlet spins coupled by J af x , the 
string order has its maximum value 0.25. With the dis- 
appearance of the singlet spins for J af x = 0, the string 
order also vanishes. In the spin- 1/2 F-AF alternating 
chain, the gap is also related to the singlet dimers coupled 
by the AF interactions, and the string order could mea- 
sure the singlet correlation Q . Hida Q pointed out that 
such a string order could distinguish the valence-bond- 
type disordered states from other disordered states. In 
this tetrameric chain, the valence-bonds, namely the sin- 
glet states in the spins coupled by Jaf ± , is characterized 
by the string order. 

The finite string order reveals the hidden Neel order 
in the gapped states of the tetrameric chain, which is re- 
garded as a feature of Haldane phase. For this tetrameric 
model, the valence-bond state could provide us a picture 
for the hidden Neel order. According to our definition of 
the string order, the exponent part of the string operator, 
which is between the spins at two boundaries, could be 
regarded as those constructed by the spin pairs coupled 
by J af 2 or Jp. Consequently, the valence bonds form 
between spins coupled by Jaf x , leading to that the spin 
pairs coupled by Jaf 2 or Jp could have Stot = 1,0,-1. 
After removing the spin pairs with Stat = 0, we can see 
a perfect Neel order of the spin pairs. 
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IV. DUALITY AND HIDDEN SYMMETRY 

In the Haldane phase of the spin-1 Heisenberg chain 
[H and spin- 1/2 F-AF alternating chain (35|, the hidden 
symmetry measured by string order was found to be a Z 2 
symmetry after a dual transformation, and the string or- 
der is transformed into a ferromagnetic order to measure 
such a Zi symmetry in the dual space. In general, the 
string orders in both x and z axes are considered to con- 
struct a Z2 x Z2 symmetry of the transformed Hamilto- 
nian. The hidden symmetry could be revealed by means 
of the dual transformation |35|,[36[ and the Haldane phase 
was characterized by a fully breaking of the Z2 x Z2 sym- 
metry 0, [32j . It is interesting that this tetrameric chain 
also shows a hidden Z2 x Z2 symmetry, which is just mea- 
sured by the string order defined above and fully breaking 
in this gapped phase. For simplicity, we start with the 
Hamiltonian rewritten in terms of Pauli matrices 

H = y~^(«/AFiO"4j-3 ' cr 4j-2 + JAF 2 <7 4j-2 ' Q~4j-1 
3 

+ JAFi^ij-l ■ &4j ~ JfC4j ' C74j + i). (22) 

After applying the standard Kramers- Wannier dual 
transformation D [37[ to the whole system, we get the 
transformed operators 

DVjGj+lD 1 = — CT j-l + l/2 cr J+l/2°"|+l + l/2; 

Da^a^D- 1 = a* +1/2 , (23) 

where c" +1 / 2 ( a — x i Vi z ) are the Pauli operators in the 
dual space. In order to present the spins more clearly, 
we merely change the labeling of the lattice sites by the 
following rule: 

#:r^i(r + 2-~). (24) 

After relabeling, we write the Pauli matrices on the 
half-odd-integer spins as r and apply the inverse of the 
Kramers- Wannier dual transformation only on r spins. 
As a result, the odd sites of both a and r spins are 
rotated by tt about the x axis. Thus, the transformed 
Hamiltonian has the form of 

H = -Jaf x /Xtftf+l + T j T j+l + a j a j+i T j T i+i) 

3 

+ ^2( J AF 2 <J%3 + Jf^j + i) + ^2(JAF 2 T%j + JfT% j + 1 ) 

3 3 

+ ^(JFOjj+iifj+i ~ ■> -u ■.■*■!>, ~L '■■ ( 25 ) 
3 

The tetrameric Heisenberg chain is transformed into a 
quantum Ashkin- Teller (AT) model in a transverse field 
which is described by Eq. P5|) . This AT model consists 
of two Ising chains coupled by four-component interac- 
tions. The transverse field parts measured by Jaf 2 an d 



Jp in this model play the role of temperature in classi- 
cal systems. They compete with Jaf x to determine the 
behavior of the system, which has been demonstrated 
in our numerical results. Similar to the quantum Ising 
model 26], a phase transition, which is referred to as a 
spontaneous breaking of the Z2 symmetry and measured 
by the ferromagnetic order, may happen because of these 
competitions. The symmetry of this quantum AT model 
can be easily read off. It is invariant under rotations of 
7r about the x axis that are applied to a spins alone or r 
spins alone. Thus, it has a Z2 x Z2 symmetry. 

After the same dual transformation, the difference of 
the four string operators defined above is revealed. The 
vanishing string orders are transformed into 

Uet^U- 1 = -® 2 fe i + 2 I +1 < (26) 

where U denotes the whole transformation. These two 
string operators are still expressed in terms of nonlocal 
operators and could not measure the spontaneous break- 
ing of the discrete symmetry. In contrast, the nonvan- 
ishing string orders become ferromagnetic correlations in 
the dual space, 

uei^u- 1 = {-iy- i+x r^ j+1 , (27) 

where the factors (— l) J ~ l and (— come from the 

rotation transformation on the odd spins. These order 
parameters measure possible spontaneous breaking of the 
Z2 symmetry of t spins. The ferromagnetic order of a 
spins corresponds to the string order defined in the x 
component. When both Jaf 2 ~ an d Jf — 0, the AT 
model becomes two ferromagnetic Ising chains coupled 
by the four-component interactions. The ferromagnetic 
correlation has the maximum 1, and the Z2 x Z 2 sym- 
metry is totally broken because of the full polarization 
along the z axis. In contrast, when JaFx = 0, the AT 
model becomes isolated rungs in a transverse field. The 
ferromagnetic correlation is 0, and the Z 2 x Z2 symmetry 
is preserved. In our model, we are considering a SU(2) 
symmetric model, thus the x component of the string or- 
der equals to that of the z component, namely O x = O z . 
It turns out that, O z could also be used to measure the Z 2 
symmetry of a spins. The nonvanishing string order O z 
indicates that the Z2 x Z2 symmetry are fully breaking 
in the tetrameric system, which is an important evidence 
for the Haldane-like type of the gapped phase. 

In this quantum AT model, the different effects of the 
couplings J af 2 and Jf could be seen more clearly. As 
shown in the transformed AT Hamiltonian Eq. (|25p . 
the odd-site spins have an antiferromagnetic interaction 
JF&2j+i T 2j+ii which has an inverse effect on the trans- 
verse field part Jf{^2q+i + r 2j+i)' an d thus counteracts 
the polarization of the odd-site spins on the x direc- 
tion. In contrast, the even-site spins have a ferromag- 
netic interaction —JAF 2 &2j T 2ji which promotes the po- 
larization on the x direction. This difference could partly 
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• spin up O spin down valence bond 



FIG. 7: (Color online) A schematic valence bond picture for 
the gapped ground state of the tetrameric Heisenberg chain, 
(a) A typical configuration of the valence bond state for the 
tetrameric chain in the absence of a magnetic field. The solid 
bullet denotes spin up, the hollow bullet denotes spin down, 
and the short line denotes the singlet bond, (b) The nonvan- 
ishing string order, (c) The vanishing string order. 



explain why the gap and string order would vanish when 
J af 2 ~^ 00 > but converge to finite values when Jp—>oo. 



VALENCE BOND GROUND STATE FOR 
THE HALDANE-LIKE PHASE 



The nature of the spin-1 Heisenberg AF model can be 
demonstrated by the AKLT model, where a Hamiltonian 
with the valence bond is constructed Q. The Haldane- 
like gapped phase of this tetrameric system implies that 
it is reasonable to expect a valence bond ground state 
that could explain the numerical results and support our 
analysis. As discussed above, the gap and string order 
are related to the singlet states of the spins coupled by 
Jaf x ■ Therefore, a valence-bond-solid (VBS) state pic- 
ture for the ground state of the tetrameric chain can be 
proposed, as shown in Fig. [7] (a). The singlet valence 
bonds, which are represented as short lines, form between 
the spins coupled by JaFi ■ The gap, which is the energy 
needed to break the bonds, should increase with J A f x ■ 
This has been confirmed by our numerical results. After 
applying a magnetic field, the magnetization plateau at 
m = 1/4 could appear when the bonds are broken, and 
the width of the plateau is thus mainly determined by 
Jaf 2 , as analyzed in Sec. III. In Figs. [7] (b) and (c), the 
exponential parts of the two topologically different string 
orders are encircled by dashed lines. The numbers under 
the lines denote the spin (5? + Sf, x ). The nonvanishing 
string order has the hidden Neel order, as shown in Fig. 
[7] (b). But in Fig. (c), the vanishing string order does 
not exhibit the hidden Neel order in this picture. 

In the absence of a magnetic field, we propose a trial 
wavefunction |f2), which is defined as a linear combina- 
tion of two functions |*) and |$), i.e. |f2)=A|*)+B|$), 
for this VBS state picture. The wavefunction Ity) is de- 



fined as 

I*) = 1> ai <8> V>a 2 <8 • • • ®VW-i ® VW£ Q1Q2 • • • e Q "- lQ " , 

(28) 

where ip ai (a, = 1,2) denote the eigenstates of S z with 
eigenvalues 1/2 and —1/2 of the spin on the site i and 
e a/3 is the antisymmetric tensor with e 12 = l/v2- These 
wavefunctions are written under periodic boundary con- 
ditions and thus L is even. This function is one of the 
ground states of the Majumdar-Ghosh model and qual- 
itatively measures the ground states of this tetrameric 
system in large Jaf ± limit. In the state \^f), the sys- 
tem has the maxima for both the gap and the string or- 
der. The singlet correlation of the spins coupled by Jaf x 
also has the maximum. In order to give a trial function 
that could depict the ground state, we introduce another 
wavefunction |$), which is defined as 



<8> 1pa 4 



aiL - 2aiL - 1 (g>ip aiL (g>ip ai X a4La \ (29) 



where A"' 3 is the symmetric tensor with A 11 = A 22 = 
and A 12 = 1/V2- |$) measures the states in both large 
Jaf 2 an d Jf limit. In this state, the system is gapless 
and the string order is vanishing. 

In the thermodynamic limit L — > oo, and |$) are 
orthogonal, namely 1 <E») = 0. The normalization of the 
trial function thus requires A 2 + B 2 = 1. We have 
performed calculations of the physical properties on this 
trial function in the thermodynamic limit. It is found 
that the trial function |f2) has vanishing local magnetic 
moment (S?) and a short-range spin correlation function. 
The string order is found to be O z string = jA 2 , which is 
finite in the state |f2) and could well fit into the numerical 
results shown in Figs. [6] (a) and (b) when the coefficient 
A is chosen as 



A 2 



x ( 



4 — Jaf 2 /Jaf! 



1. 



3^4 - 2J AF JJ AFl + (Jaf 2 /Ja Fi ) 2 3 

4 + Jf/J AFi i } (30) 

3^4 + 2J f /Ja Fi + (Jf/Ja Fi ) 2 3 y 



The numerical results of the gap as well as this valence- 
bond-state picture indicate that the gap is induced by the 
valence bond, namely the singlet state between the spins 
coupled by Jaf x - The string order could well describe 
the gap and the singlet correlation that measures the 
singlet state of the spins coupled by Jaf x ■ This above 
picture shows the crucial role of the valence bond in the 
formation of the Haldane-like gap and hidden symmetry. 



VI. SUMMARY AND DISCUSSION 

By means of the DMRG method, we have studied the 
magnetic properties, spin gap, spin-spin correlation func- 
tion, and string order of the spin- 1/2 tetrameric HAFC 
with alternating couplings AF 1 -AF 2 -AF 1 -F. Two mag- 
netization plateaux at m = and 1/4 are found, which 
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satisfy the OYA condition for such a spin chain with a 
period of 4. For the m = 1/4 plateau states, an ap- 
proximate wave function is proposed, which fits into the 
numerical results of the local magnetic moment and spin 
correlation function perfectly. The magnetization process 
of the two regions between the plateaux is fitted well by 
Eq. ©. 

Besides the m = 1/4 plateau, the gapped ground state 
in the absence of the magnetic field is also investigated. 
The system is found to have a gap from the singlet ground 
state to the triplet excited states. The gap changes mono- 
tonically with the couplings and does not show any level- 
crossing. It decreases with Jp and Jaf 2 , but increases 
with Jaf x - The ground state energy has no nonanalyt- 
icity in the parameter space, which implies an absence 
of the quantum phase transition. The asymptotic behav- 
iors of the gap have been calculated, which have the same 
features as those in the isolated tetramer systems. Com- 
bining the numerical results with the nonlinear a model, 
it is shown that a gap would be generated by an arbi- 
trary small Jafi an d thus this system is always gapped 
if J AFl ^0. 

The spin-spin correlation function is uncovered to de- 
cay exponentially in the gapped states, but the corre- 
lation length and the behaviors of the correlation func- 
tion show dramatic changes for different couplings, which 
are observed in the variations of the peaks of the static 
structure factor S(q), reflecting a complex competition 
of the interactions. When J af 2 or Jf is small compared 
with Jafi j S(q) is dominated by short-range correlations. 
It would exhibit three maxima at q = tt/3, tt and 57r/3 
when Jp is small. When Jaf 2 is small, the maximum at 
q = tt becomes a valley and two small maxima appear 
near q = Sir/ '5 and 7ir/5. In contrast, when both Jp 
and Jaf 2 are large enough, S(q) exhibits four peaks at 
q = 7r/4, 37r/4, 5tt/4 and 77r/4, which is due to the in- 
crease of the correlation length and the translation sym- 
metry of the ground states. 

In order to investigate the hidden symmetry of the 
gapped ground states, we proposed a string order and 
performed calculations to detect the hidden symmetry. 
The string order is found to be nonvanishing in the 
gapped phase, indicating a hidden symmetry. More im- 
portantly, the behaviors of the string order have the same 
tendencies as that of the gap. The opening of the gap 
and its variations could be well described by this nonlo- 
cal order. The gap is related to the singlet state of the 
spins coupled by J 'aFi &ud increases with this interaction. 
Thus, the gap could be regarded to be determined by the 



singlet correlation. The nonlocal string order could also 
measure this local singlet correlation. The singlet state 
of the spins coupled by J AFi provides a picture to under- 
stand the hidden Neel order reflected by the finite string 
order. Different from the conventional hidden Neel order, 
the order in the present system is not of the singlet spins 
but of the spin pairs coupled by Jaf 2 or Jf. 

The hidden symmetry measured by the string or- 
der is unveiled by a dual transformation. The present 
tetrameric Hamiltonian is transformed into a quantum 
AT model, which has a discrete Z^ x Zi symmetry. The 
string order becomes a ferromagnetic order that is proper 
to measure this discrete symmetry. The nonvanishing 
string order indicates that this Zi x Z2 symmetry is fully 
breaking in the gapped phase, which is the characteris- 
tic of Haldane phase. Thus, this tetrameric chain is in a 
Haldane-like gapped phase. 

A valence-bond-solid state picture and a trial function 
for the gapped states are proposed, most properties of 
which could be explained on the basis of this picture. 
The numerical results as well as the VBS picture indicate 
that the gap is induced by the valence bond of the spins 
coupled by Jafi ■ The valence bond is measured by the 
singlet correlation. The string order could well describe 
not only the gap, but also the singlet correlation. The 
critical role of the valence bond on the formation of the 
Haldane-like gap and hidden symmetry is unveiled. 

Finally, we would like to state that unlike the spin- 1/2 
F-AF alternating HAFC, the present spin- 1/2 tetrameric 
HAFC could not reduce to an integer-spin HAFC. There- 
fore, our findings extend the substance of Haldane's con- 
jecture that was originally proposed for HAFCs with inte- 
ger spin. We expect that our investigations would deepen 
further understanding on the physical properties of low 
dimensional quantum magnetism. 
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